AD-A233  938  G> 

NASA  AVSCOM 

Technical  Memorandum  102422  Technical  Report  89-C-013 


A  Planar  Reacting  Shear  Layer  System 
for  the  Study  of  Fluid  Dynamics- 
Combustion  Interaction 

C.J.  Marek  and  C.T.  Chang 
National  Aeronautics  and  Space  Administration 
Lewis  Research  Center 
Cleveland,  Ohio 


B.  Ghorashi  and  C.C.  Wey 
Cleveland  State  University 
Cleveland,  Ohio 

C.  Wey 

Sverdrup  Technology,  Inc. 

NASA  Lewis  Research  Center  Group 
Cleveland,  Ohio 


FAv  -pr 
8  fl  1  |  t 

X  s'. 

Mk  EXECTE 

MAR  2  7  1991  i  j$ 

W  <  h  i’j 


E.J.  Mularz 

Space  Propulsion  Directorate 

U.S.  Army  Aviation  Research  and  Technology  Activity— AVSCOM 
Lewis  Research  Center 
Cleveland,  Ohio 


Prepared  for  the 

23rd  International  Symposium  on  Combustion 
sponsored  by  The  Combustion  Institute 
University  of  Orleans,  France,  July  22-27,  1990 


SYSTEMS  COMMAND 
AVIATION  R4T  ACTIVITY 

'  T  A  i 


91  3  22  080 


A  PLANAR  REACTING  SHEAR  LAYER  SYSTEM  FOR  THE  STUDY  OF  FLUID 


DYNAMICS-COMBUSTION  INTERACTIONS 

C.J.  Marek  and  C.T.  Chang 
National  Aeronautics  and  Space  Administration 
Lewi s  Research  Center 
Clevel and,  Ohio  441 35 

B.  Ghorashi  and  C.C.  Wey 
Cleveland  State  University 
Cleveland,  Ohio  441 15 

C.  Wey 

Sverdrup  Technology^  Inr- 
NASA  Lewis  Research  Center  Group 
Cleveland,  Ohio  44135 

and 

E.J.  Mularz 

Propulsion  Directorate 
U.S.  Army  Aviation  Research  and  Technology  Activity  -  AVSCOM 

Lewi s  Research  Center 
Cleveland,  Ohio  44135 


/ 


ABSTRACT 


A  versatile  planar  reacting  shear  layer  (PRSL)  facility  is  constructed  at 
NASA  Lewis  Research  Center.  This  paper  describes  the  research  objectives,  as 
well  as  design,  instrumentations  and  the  operational  procedures  developed  for 
the  system.  The  fundamental  governing  equations  and  the  type  of  quantitative 
information  that  are  needed  from  experiments  are  described.  Additionally,  a 
review  of  earlier  work  is  presented.  Whenever  appropriate,  comparisons  are 
made  with  similar  systems  in  other  facilities  and  the  main  differences  are 
described.  Finally,  the  nonintrusive  measurement  techniques;  PLIF,  PMS,  LDV 
and  Schlieren  photography,  and  the  type  of  experiments  that  are  planned  are 


described . 


INTRODUCTION 


Several  recent  studies  and  reviews  of  current  research  in  turbulent 
reacting  f:ows  have  concluded  a  clear  need  for  additional  work  in  plane  free 
shear  layer  systems. ^  Although  a  repertory  of  literature  is  available  on  round 
free  jets,  relatively  few  works  are  reported  on  plane  free  shear  layers.  At 
the  same  time,  there  is  an  imperative  need  for  further  understanding  of  the 
plane  mixing  layer  in  order  to  develop  relationships  between  the  existing 
computational  fluid  dynamics  codes  and  the  experimental  results. 

For  the  above  reasons  and  in  order  to  provide  an  understanding  of  fluid 
aynami cs-cembustion  interaction  for  space  and  aeronautics  projects,  a  very 
substantial  and  versatile  planar  reacting  shear  layer  f PRS L )  system  is  being 
designed  and  built  at  NASA  Lewis  Research  Center. 

The  importance  of  a  planar  reacting  shear  layer  (TflSL)  facility  to  NASA 
Lewis  Research  Center  can  be  described  in  terms  of  space  and  aeronautics 
applications.  For  space  applications,  there  is  a  need  for  development  of 
methods  to  predict  low  thrust  propulsion  systems  in  terms  of  specific  impulse 
and  wall  heat  transfer  as  two  examples.  The  design  methodology  of  large 
rocket  motors  poorly  predicts  low  thrust  systems.  For  instance  the  hydrogen 
oxygen  mixing  layer  of  thruster  can  he  simulated  with  a  PRSL  system  (Fig.  1) 

For  aeronautics  applications  (i.e.,  high-speed  propulsion),  a  PRSL  system 
would  mate  it  possible  to  examine  the  effect  of  acoustic  waves  and  the  effect 
of  combustion  on  mixing  (e.g.,  gas  temperature  rise)  and  engine  nozzle 
^M^,„.dMce  as  well  as  the  effect  of  mixing  on  combustion  intensity  (Fig.  2). 
Additional  benefits  of  this  facility  include  the  capability  to  conduct 
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extensive  high  temperature  diagnostics  for  heat  transfer  research  and  new 
materials  evaluation.  The  main  objectives  are  (1)  to  investigate  the 
mechanisms  that  control  the  mixing  of  fuel  and  air,  centerline  variations  in 
mixture  fraction,  mixture  fraction  fluctuation  and  intensity,  (2)  to  predict 
the  shear  layer  spreading  rate  and  the  effects  of  the  initial  conditions  on  jet 
development,  (3)  to  determine  the  effect  of  large  scale  structures  on  mixing  as 
well  as  the  effect  of  excitation  on  the  structures,  (4)  to  identify  the  effect 
of  heat  release  on  turbulence  flowfield,  and  (5)  to  more  closely  examine 
variable  density  flows  (i.e.,  density  difference  between  fuel  and  coflowing 
air). 


Products  of  Experiments 

Specifically,  the  products  of  experiments  would  be:  (1)  data  for 
validation  of  emerging  computational  fluid  dynamics  codes  (CFD), 

(2)  measurement  of  interaction  between  flow  acoustics  and  combustion  processes 
and  (3)  insight  in  exploring  mixing  enhancement  techniques.  There  will  be  an 
emphasis  placed  on  pressure  wave-shear  layer  interaction  in  order  to  determine 
the  coupling  between  fluid  dynamics  and  combustion,  to  validate  CFD  codes  and 
to  accurately  simulate  the  physics  and  chemistry  of  high  speed  chemical 
reacting  flow.  This  requires  an  unsteady  two-dimensional  flow  code  with 
pressure  oscillation-fluid  interactions,  combustion-pressure  interactions  and 
finite  chemical  kinetics.  If  ’s  determined  that  in  order  to  realize  the  above 
goals,  the  following  information  should  be  obtained:  (1)  a  correlation  between 
pressure  oc-r ;  1  1  ations  and  larae  pddv  fnrrnaHon.  (2)  tiara  al  rnui,  neai  release 
in  order  to  correlate  the  combustion  rates  to  eddy  growth  rates  and  pressure 
changes,  (3'  measurements  in  high-speed  compressible  flows,  and  (4)  high 
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pressure  data  by  using  a  downstream  nozzle  to  determine  the  effects  of 
oscillating  pressure  field  on  shear  layers. 

Background  and  Early  Studies 

As  stated  earlier,  the  available  literature  on  plane  free  shear  layers  is 
comparatively  limited  and  what  is  available  is,  in  most  part,  restricted  to 
simple  turbulen4:  flows  with  no  chemical  reactions.  Furthermore,  virtually  all 
of  these  studies  are  limited  to  simple  geometries  without  the  complication  of 
radiative  heat  transfer  and  multiple  phases.  It  is  further  understood  that 
the  above  complications  can  and  will  have  profound  and  drastic  effects  on  the 
final  outcome  and  therefore  limit  the  degree  of  suitability  of  the  results  of 
those  simpler  studies. 

Among  the  works  on  turbulent  reaction  in  plane  free  shear  layers,  under 
the  fast  reaction  limit  studies,'  Mungal,  et  al.2  studied  a  shear  layer  in  a 
duct  with  dilute  hydrogen/f 1 uor  i ne  in  nitrogen  or  helium  under  the  conditions 
of  Re*  =  4 x 1 0 ^  and  132  K  maximum  temperature  rise.  Koochesfahani  and 
Dimotakis2  also  investigated  a  shear  layer  in  a  duct  with  dilute  acid/base 
reaction  in  water.  They  used  H2SO4  and  NaOH  with  a  time-resolved  LIF 
visualization  along  a  line  crossing  the  flow.  Batt4  investigated  a  wall  jet  in 
still  air  via  dilute  nitrogen  tetroxide  dissociation.  Wallace2  studied  shear 
layer  in  a  duct  with  dilute  nitric  oxide/ozone  in  helium,  nitrogen  or  argon. 
Also  many  studies2  are  reported  on  means  to  control  shear  layers.  These  are 
mostly  controlled  CaC  ;  Ld i  u"  i.echn  I  goes  applied  Lo  noi  u  e  a  c  t  i  m  g  sneai  layers 
with  very  promising  results. 


4 


SYSTEMS  DESCRIPTION 


The  planar  reacting  shear  layer  features  a  continuous  flew  under  high 
subsonic  velocity  conditions  and  high  heat  release.  Figure  3  depicts  the 
overall  scnematic  of  the  system.  In  order  to  distinguish  the  different 
features  of  this  facility,  the  description  that  follows  is  divided  into  three 
sections,  namely:  (1)  the  PRSL  Test  Section,  (2)  cne  Laser  Diagnostics  and 
Pressure  Waves  Instrumentation  System,  and  (3)  the  Image  Acquisition  and 
Processing  System. 


PRSL  Test  Section 

The  cross  section  of  the  test  section  is  10.2  by  20.4  cm  at  the  location 
of  the  knife  edge  (Fig.  4).  The  height  (20.4  cm)  is  large  in  order  to  keep  the 
shear  layer  from  being  influenced  by  the  boundary  layers  on  the  wall.  The 
interaction  between  N£,  H£,  and  heated  air  (867  K)  begins  at  the  knife  edge 
This  section  is  preceded  by  screens  and  a  honeycomb  section  for  fuel  and  air. 
The  reacting  shear  layer  discharges  into  a  larger  chamber  and  passes  over 
several  water-cooled  back  pressure  tubes  (Fig.  5).  These  tubes  are  removable 
in  part  or  in  their  entirety.  Microphone  pressure  taps  are  placed  along  the 
upper  and  lower  sides  of  the  test  section  in  order  to  investigate  the 
fluctuating  pressure  field.  For  nonintrusive  measurements,  laser  sapphire 
windows  are  constructed  on  the  sides.  Sapphire  is  a  superior  window  material 
because  of  its  extreme  surface  hardness,  high  electric  constant,  high  thermal 
conductivity  and  resistance  to  most  chemical  acids  and  alkalis.  Because  of 
their  structuia!  strength,  sapphire  windows  can  be  made  much  thinner  than  other 
materials.  Chemically,  sapphire  is  a  single  crystal  a  1 umi num  ox i de  ( A 1 2 O3 ) 
and  exhibits  high  internal  light  transmittance  from  150  to  6000  nm.  The  upper 
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and  lower  walls  are  moveable  la  order  to  be  able  to  impose  a  desired  axial 
pressure  gradient. 


Laser  Diagnostics 

The  most  important  tasks  are  to  identify  the  instantaneous  structure  of 
the  flowfield,  the  location  of  the  burning  zones  and  the  effect  of  finite-rate 
cnemistry.  It  should  be  noted  that  although  planar  shear  layer  is  the  simplest 
shear  layer  geometry  and  is  well  suited  for  modelling  purposes  (through 
numerical  computations),  it  is  experimentally  difficult  to  implement. 

The  following  factors  are  of  paramount  importance  in  choosing  diagnostic 
systems : 

I.  Instantaneous  structure  of  mixing  zones.  -  Mixing  begins  to  occur  as 
the  fuel  and  oxidizer  (air)  come  together  at  the  knife  edge.  This  process 
cannot  be  inferred  from  time-averaged  concentration  measurements  since  those 
techniques  can  only  depict  an  overall  picture  of  the  mixing  process.  A  more 
detailed  knowledge  of  the  mixing  process  and  the  flowfield  can  be  obtained  by 
measurements  of  the  instantaneous  concentration  field.  This  understanding  can 
ce  further  enhanced  by  a  juxtaposition  of  the  instantaneous  velocity  and 
pressure  fields. 

II.  Location  of  burning  zones.  -  A  detailed  study  of  burning  zones  is 
essential  in  order  to  construct  a  relationship  between  the  temperature  fields 
and  the  local  velocity  fields  Also,  it  is  equally  important  to  investigate 
she  differences  between  burning  under  supersonic  and  subsonic  conditions. 
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III.  Effect  of  finite-rate  chemistry.  -  it  is  anticipated  that  the 


Damkohler  number  (mixing  time/chemical  reaction  time)  would  be  small  in  certain 
regions  of  the  flow. 2  f h i s  finite-rate  chemistry  may  play  a  significant  role 
in  supersonic  combustion.  The  required  information  on  the  effect  of  finite- 
rate  chemistry  may  be  obtained  by  measurement  of  intermediate  and  product 
species  as  well  as  temperature. 

V.  Two-component  LDV  measurement.  -  A  two-component  (four-beam)  LDV  system 
is  used  to  make  point-measurement  of  the  time-evolved  velocity  fluctuations 
inside  the  flowfield.  In  order  to  obtain  a  hi gh- i ntens i ty  scattered  signal 
for  high  signal  to  noise  ratio,  a  3.75x  beam  expander  and  forward  scattering 
are  used.  However,  this  reduces  the  possibility  of  particle  detection  and 
decreases  the  data  rate. 

Conventional  Instrumentations  and  Advanced  Diagnostics 

LDV,  pressure  probe  and  Schlieren  photography  are  among  the  conventional 
instrumentations  that  will  be  used  for  measurements  that  are  described  above. 
Among  the  advanced  diagnostic  techniques,  planar  laser  induced  fluorescence 
(PLIF)  will  be  used  for  measurements  of  the  instantaneous  concentration, 
temperature,  velocity  and  pressure  fields.  Also,  planar  Mie  scattering  (PMS) 
will  be  used  to  obtain  the  instantaneous  physical  image.  Figure  6  is  a  chart 
showing  the  schematics  of  different  PLIF  systems  and  their  applications  for 
temperature  and  concentration  measurements.  Other  schematics  on  the  chart  are 
velocity  and  pressure  imaging  and  PMS  systems.  Figure  7  is  a  schematic  of  the 
location  of  the  laser  sheet  in  the  facility.  The  motivations  for  performing 
LDV  measurements  are  to  obtain  time-averaged  reacting  flow  data  of  the  shear 
layer  (e.g.,  mean  and  root  mean  square  U  and  V  throughout  the  flowfield) 
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as  a  benchmark  for  comparison  with  computation  models,  and  secondly,  to 
examine  the  effect  of  oscillating  pressure  field  on  shear  layer  through  the 
f'oss-correl ations  of  instantaneous  pressure  and  velocity.  Table  1.5  As 
mentioned  earlier,  the  objectives  are  to  obtain  time-dependent  data  of  mean 
and  turpulence  quantities  as  well  as  images  for  comparison  with  predictions. 

In  order  to  increase  the  data  rate  and  obtain  high-density  data,  a 
reactive-type  seeder  is  "sed  to  generate  seed  density  of  about  one,  0.3  urn 
diameter  T i O2 .  particle  per  cubic  millimeter.  These  particles  are  formed  from 
the  reaction  of  small  droplets  of  titanium-tetrachloride  and  steam.  In  case 
tne  scatter  cross  section  of  titanium  dioxide  becomes  too  small  within  the  higi 
temoerature  gas  trimethyl-aluminum  can  be  used  as  a  direct  substitute  for 
titanium-tetrachloride.  The  formed  alumina  particles  are  quite  stable  in  high 
temperature  conditions.  This  chemical  particle  generation  method  yields  a  more 
e.en  distribution  and  higher  density  than  the  dry  fluidized  bed  seeders. 

Image  Acquisition  and  Processing  System 

A  scnematic  diagram  of  the  imaging  system  is  illustrated  in  Fig.  8.  This 
system  is  comprised  of  a  host  computer,  an  image  processing  subsystem,  image 
acquisition  hardware  and  data  analysis  and  reduction  components.  The  cameras 
are  of  the  gated  intensified  solid-state  CCD  monochrome  video  camera  and  can  be 
t-iggered  externally.  PLIF  requires  one  camera  for  concentration  measurements 
.i'b  two  sameras  for  temperature  measurements.  Velocity  measurements  require 
'wo  came-as  for  one-dimensional  and  four  cameras  for  two-dimensional 
measurements.  Also,  since  the  fluorescence  intensity  is  too  low  to  be 
measured  directly,  a  UV-enhanced  microchannel  plate  (MCP)  image  intensifier  is 
internally  bonded  to  the  CCD  imaging  array  via  an  optical  fiber  bundle.  This 
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system  offers  488(V)  by  754(H)  pixel  frame,  30  f rames / sec .  An  image  subsystem 
that  controls  a  video  tape  recorder  through  an  animation  controller  is  used  to 
acquire  images.  For  image  processing  and  analysis,  single  video  frames  wil; 
oe  acquired  from  the  previously  recorded  tape  and  will  be  fed  into  the  image 
subsystem. 

The  host  computer  is  a  Concurrent  5600  (Fig.  9).  The  operating  system  :s 
a  real-time  enhanced  UNIX  system  (RTU  4.0)  that  is  compatible  with  AT&T  System 
V  and  Berkely  4.2  BSD. 


THEORY  AND  COMPUTATIONS 


The  Large  Scale  Coherent  Structure  and  Initial  Condition  Effects 

The  PRSl  is  a  two-dimensional  planar  variable-density  flow  system.  It 
should  be  noted  that  the  term  two-dimensional  is  soley  descriptive  of  the 
system  geometry  as  any  turbulent  flowfield  is  intrinsically  three-dimensional 
in  nature.  As  such,  the  large-scale  structures  play  an  important  role  in  the 
mixing  process.  In  round  jets,  density  variations  occur  in  regions  downstream 
of  the  flow  and  decrease  with  increase  in  the  axial  distance.  On  the  other 
hand,  two-dimensional  shear  layers  offer  the  advantage  of  providing  density 
differences  farther  downstream  (i.e.,  a  longer  mixing  region).  Therefore,  any 
modeling  attempt  should  incorporate  the  effects  of  large  scale  structures  on 
the  mixing  process.  Zaman  and  Hussain^  have  summarized  the  results  of  their 
studies  on  natural  large-scale  structures  in  the  axisymmetric  mixing  layer. 
They  have  also  made  comparisons  between  the  natural  structure  and  the 
structure  induced  by  controlled  excitation.  Among  their  findings,  they 
report  that  the  natural  structures  were  independent  of  Rep  over  a  range  of 
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5.5<10*  y  Reg,  y  3  <  1 0 D .  Tne  educed  structures  for  the  natural  and  excited  cases 
.re  .  y  $ ; m* ' a-  except  at  nigncr  amplitudes  where  excitation  can  noticeab’y 
a*  f-e  ;:rcr  er  -or  tne  analytical  f'  emulation  on  the  treatment  of  the 
:  "e-  ert  a-i  --_--.ro.-ort  f  lcwf  ie  ids ,  the  reader  may  consult  Ret.  6. 

.Is mar  arc  •  . ; s a ' n  a’so  report  a  rila  dependence  on  the  initial  boundary 
’  i  y  e  ’■  to-d !  ts  c-  wi  tn  lan'nar  boundary  1  ^>er  producing  more  disorganized 

re  s .  :  s  ’  atte-  -’nding  is  also  supported  by  Clark.  ;  Hussain  ana 

Us "a-  1  snow  :-e  contours  of  coherent  spanwise  vorticity  at  different  stages 
•'*  :■-$  :  t  -  j  c  t  u  r  o  evolution  in  a  p'ane  mixing  layer  (Fig.  10);  the  last  station 
. :  '  e see- d s  :.  a  Reynolds  number  of  Rex  ~  2x10^.  The  .ontours  are  plotted  as 
a  -  met  i  3-  ~,f  tie  transverse  distance  and  the  time.  The  authors  note  that 
w-r'e  t-e  st-wCt ure  shape  cnanges  between  x  =  150  and  300  6e,  the  change  i,i 
t-e  rona!rar;:rr,a’  vorticity  is  not  significant.  They  also  observe  that 
'  a-  •:  e-  c  /c-t  i :  •  ty-tontcur  details  suggest  achievement  of  an 

eg  s -  1  ’  o- 1  rs  state  whizn  persists  for  the  enfre  length  of  plane  mixing 
: a  ye r .  1 


Effect  of  Heat  Release  on  Tu-bulent  Flowfield 


rarsh:nid  recently  extended  a  fully  second-order  closure  model  for 
■easting  turbulent  flows  in  order  to  include  the  effects  of  heat  release  on 
turbulent  f’owfield.  He  analyzed  the  exact  pressure  equation  of  a  turbulent 
-enpremixed  flame.  His  completely  modelled  scalar  d;ssipation  equation  is: 
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F a r s h c h i  made  comparisons  with  limited  experimental  results  but  noted  that 
furtner  experiments  are  needed  to  support  or  to  improve  his  assumptions.  He 
specifically  suggested  comprehensive  measurements  of  quantities  such  as  mixture 
fraction,  ocnsity  variances  and  scalar-velocity  correlations.  These  fall 
exactly  within  the  scopes  and  objectives  that  are  planned  for  the  PRSl  system. 


The  modeled  Reynolds  stress  equation1®  for  the  case  of  heat  release  with 
the  assumption  of  local  isotropy  at  high  Reynolds  numbers  is  as  follows: 


■ , ( p  u,  Uj'uvJ  =  -(o  udut  a] +  P  a1  Uj 


-  3  P  c  Sjj 


Cr  p  -  uVu"  9  (  u  •'  u 
b  e  I  m  mV  i  i 


+  -  i  p 1  u 9  .  P  x-  p 1  u 3.P)  +  0-  •  x-  0  -  • 

P  i  ]  H  J  i  /  ij  vij 
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where  represents  the  correlation  of  the  velocity  f 1 ucturat ions  with  the 

hyarcdynami cs  part  of  the  pressure  fluctuation  gradient  and  represents 

the  heat  release  part  of  the  pressure-strai n  rate  correlation . ^ 

Measurment  of  the  kinetic  energy  term  (k)  is  of  particular  importance 
since  it  :dent’fies  the  level  of  turbulence.  As  such,  the  following  quantities 
are  of  primary  ’importance: 

u"  ,  U'.'ut  k.e  ,  p'u" 

CONCLUDING  REMARKS 


A  detailed  description  of  a  versatile  planar  reacting  shear  layer  system, 
jnder  construction  at  NASA  Lewis  Research  Center,  is  provided.  The  motivations 
arc  cojecti/es  for  the  design  and  the  type  of  quantitative  information  that 
are  needed  from  experiments  are  described.  State-of-the-art  diagnostic  systems 
incorporated  in  the  PRSL  system  and  their  interrelations  are  discussed, 
finally,  a  synopsis  of  the  theory  and  the  required  computations  is  provided 
which  includes  the  effect  cf  large  scale  structures,  influence  of  the  initial 
Dcundary  layer  conditions,  effect  of  excitation  amplitude  on  the  structures 
and  the  effect  of  heat  release  on  the  turbulent  flowfield.  The  PRSL  system 
would  enable  researchers  to  develop  relationships  between  the  existing  CFD 
codes  and  the  experimental  results  to  further  explore  fluid  dynamics/combustion 
interaction  for  space  and  aeronautics  projects. 
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NOMENCLATURE 


f  mixture  fraction,  f  =  f  +  f" 

k  kinetic  energy  of  turbulence  (1/2  UiUp 

P  Favre  probability  density  function 

p  fluctuating  pressure 

Re  Reynolds  number 

U,V  mean  velocity  components 

u,v,w  fluctuating  velocity  components 

x  position  vector 

x,y,z  streamwise,  transverse  and  spanwise  coordinate 

c  rate  of  tii ssipation 

e -  scalar  variance  dissipation  rate 

v  kinematic  viscosity 

o  fluid  dens i ty 

(  )  conventional  averaging 

(  ),(  )"  density  weighted  averaging 

ABBREVIATIONS 

CFD  computational  fluid  dynamics 

LIE  laser  induced  fluorescence 

LDV  laser  doppler  velocimeter 

MCP  microchannel  plate 

PDA  photodiode  array 

PLIF  planar  laser  induced  fluorescence 

PMS  planar  mie  scattering 

PRSL  planar  reacting  shear  layer 
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1AIUE  1.  -  DESIRABLE  QUANTITIES  FOR  MEASUREMENT1' 
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